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I/"") \ Abstract 

The static three-quark potential in arbitrary configuration of quarks is calculated 
analytically. It is shown to be in a full agreement with the precise numerical simulations 
in lattice QCD. The results of the work have important application in nuclear physics, 

\f} , as they allow to perform accurate analytic calculations of spectra of the baryons. 
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^ ; 1 Introduction 

o ; 

Precise knowledge of the nucleon spectroscopy is very important for the deep understanding 
of the nuclear processes. The nonperturbative static potential is a key quantity in calculations 
of the baryon spectra, and this is why its investigation is one of the important problems in 
nuclear physics. Besides, the study of the static potential provides an insight into the two 
fundamental properties of strong interactions — confinement of the color charges and scale 
of the fluctuations of the confining gluonic fields. 

In the present paper static potential in baryon is calculated using the method of the field 
correlators (MFC) H, which is based directly on QCD and operates with the vacuum averages 
of the correlators of the gluonic field strengths. 

It is suggested in MFC that the gluonic fields play a twofold role in QCD. First, gluons 
propagate dynamically in the vacuum, and this process at small distances can be described 
by the pertubation theory. In particular, the interaction of the quarks due to the one-gluon- 
exchange leads to the color-coulomb interaction potential. Second, gluons in the vacuum form 
the nonperturbative condensate, which constitutes the background, where the perturbative 
gluons propagate. 

Two-point (or bilocal) correlators of gluonic fields are parameterized in MFC by the scalar 
formfactors ||. This parameterization gives rise to the area law for the Wilson loop at 
sufficiently large separations between color sources, which ensures confinement. The slope a 
of the linear static potential between quark and antiquark, corresponding to the Wilson loop, is 
the main parameter of MFC. Using this the only parameter one can obtain in MFC the spectra 
of the mesons within the 5% accuracy ||. The value of a is determined phenomenologically 
from the slope of the Regge trajectories of mesons, and corresponds to the confinement radius. 

The scalar formfactors of the background fields fall off exponentially, which reflects the 
stochastic nature of the confining fluctuations of the gluonic background field. The correlation 
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length of the background fields T g is sometimes considered as the second parameter of the 
MFC, although it can generally be expressed through the o [J12||. Its value is small in compar- 



ison with the radius of confinement and negligible for mesons. The situation for baryons is 
different because of the correlations of the background fields near the string junction, which 
lead according to MFC to the decreasing slope of the potential at small and intermediate 
quark separations. There is another T 9 -induced effect, which is considered in the present 
work, namely the potential difference between the configurations with the same length of the 
string, but different locations of quarks. 

The appearence of recent accurate numerical calculations of the static three-quark poten- 
tial in lattice QCD in quenched approximation |5|, || allow to test analytic results of MFC. 
We perform this important test in the paper and find a complete agreement between MFC 
and lattice results. 

Before to calculate the static potential, in the next section of the paper we consider 
the Wilson loop of the baryon, which in the case of the static quarks is the baryon Green 
function. We express vacuum average of the Wilson loop through the vacuum averages of 
bilocal correlators of the gluon field strength tensor and parameterize them according to MFC. 
A brief discussion of the gluon condensate value obtained from the ITEP sum rules and its 
relation with the MFC is also given. In third section the procedure of calculation of the 
potential is described, and its behavior is analized in comparison with the lattice data. In 
the concluding section a summary of the results obtained in the paper is given, and their 
perspective applications are discussed. 



2 The baryon Wilson loop and its average in MFC 

It is known that the gauge-invariant state of the baryon with quarks at points x, y, z reads 
as 

* B (x,y,z,Y) = e^^{xXCi)^{yXC2)^{zXC 3 )^{x)g^{y)q^(z^ (1) 
where a, j3, ... are the color indexes of the fundamental representation of SU(3), 

^(x,y,C) = (Pe W ig [ A„dz^ (2) 

Jo 

is a parallel transporter along the contour C, which connects points x, y (P in (0) means 
the ordering of color matrices along the trajectory of integration); q and A are quark and 
gluon field operators, and Y is the point of the string junction. The baryon Green function 
is defined as a vacuum average 

G b (X,X) = (V+(X)*b(X)), (3) 

where X = x, y, z, Y. If quarks are static, the Green function becomes the Wilson loop Wb, 
which reads as 

Wb = (^ea^^^Ct)^^)^^ . (4) 

The contours C\, C^, C3 and direction of the integration are shown in Fig. 1. The 
trajectory of the string junction, which we denote C , is shown in the Fig. 1 by the dotted 
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Figure 1: Baryon Wilson loop 



line. Note that the definition (^) of the baryon Wilson loop is generally accepted. The static 
potential in baryon is related with the Wilson loop as follows, 

y B = - T lim ilnW Bl (5) 

where T is the time extension of the Wilson loop. 

We are now to express the Wilson loop (f|) through the correlators of the gluon field 
strength tensor. For this sake rewrite it as 



Wb = ( «^e a ' W C^i)<(C 2 )€(C 3 ) ) , (6) 



1 

where C{ = Ci UCV The expression (|(J) follows from (^) after the substitution the latter with 
the equation 

Note that the last relation is valid in the generalized coordinate gauge, in which $^,(Co) = S^,, 
but the expression (|6]) is a scalar of the color SU(3) and therefore does not depend on the 
choice of the gauge. Since the integration in is performed along three closed contours, one 
may use for these contours of the baryon Wilson loop the nonabelian Stokes theorem 0. As 
a result one gets 

1 
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W B = ( ^e^'^'W^W^W'^Ss) ) , (8) 



where Si are the minimal surfaces of the contours C«, and 



W«,(S) = VeMw I da, u (x)F^(x)<S>(x,x )) a a , (9) 
Js 
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does not depend on the choice of the point x , if the latter places at the surface S Expand 
now the exponent in (^), 

W%(S) = 5«, +igj s da F« -£ J J da(l) da(2) Ff(l)F*(2) + (10) 

where we have explicitely written only the color indexes and denoted by the dots the higher 
powers of the expansion. Since the vacuum average of the strength tensor equals to zero, the 
Wilson loop will read as follows, 

-Zfff I da{l)da{2) 1^(1)^(2) + ...Y (11) 
J Si Js 2 J 

where the dots denote the double integrals over surfaces S2S2, S3S3, Si S3, and S2S3, anal- 
ogous to the written ones, as well as the higher order correlators. Taking into account that 
eafrf? & F?F£ = 2trFF, e a ^e a '^F^,F^ = -tiFF, one gets from (0) 

W B = exp{-^ I I da(l)da(2) / 9 jtrF(l)F(2)\ + 

^ j = l J Si J Si \ / 

+ f da(l)da(2)/ 9 jtrF(l)F(2)\\. (12) 

i<j J Si J Sj \ /j 

Here and in what follows we disregard the higher order correlators, assuming that their 
contribution to Wilson loop is small or can be reduced to the contribution of the bilocal 
ones. This assumption leads to the area law for the Wilson loop and is shown in |7j to be 
motivated by the property of the Casimir scaling of the potentials of static sources in various 
representations of SU(3). 

The second motivation of the cancellation of the contributions of higher order correlators 
to the Wilson loop follows from the stochastic properties of the QCD vacuum. Note that 
the bilocal correlators written in (0) are the first term of the cluster expansion used in the 
theory of fluctuations [|] for the description of the correlated stochastic processes. According 
to the fundamental identity of the cluster expansion, higher order correlators enter into fli"2D as 
cumulants, or connected correlators (one can find an accurate definition of cumulants as well 
as the proof of the fundamental identity in 0), which rapidly fall off with the increasing order 
for almost independent points. That's why the bilocal approximation taking into account only 
two-point correlators is motivated by the theory of fluctuations at distances greater than the 
correlation length of the gluon fields. 

The bilocal correlators are parameterized in MFC in the case of SU(N C ) using two 
scalar formfactors D and D\ as follows, 



^\F/j,iui (x)$(;E, X )F^ 2U2 (x )$(x , x)) ifi i±iiiq^vxV2 ^ PiV2^ H2V1) F) {z} -\- 
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where z = x — x'. The formfactors of the background fields read as 



D(z) = D(0) exp ( -H j , Dl (z) = D 1 (0) exp ^-H ) . (14) 



The exponential behavior of the formfactors is confirmed by the lattice simulations [|10| , |Tlf at 
distances z > 0.2 fm. The values of the parameters obtained are T g = 0.12 -h 0.2 fm [[10], |TT 
and Di(0)/D(0) « 1/3 0. 

In the case of static quark and antiquark the relation analogous to ( |T2"D reads as 



exp { ~ JJ s da(l)da(2) ^trF(l)F(2)^ (15) 



and at large quark separations R^> T g gives rise to the area law for the Wilson loop, 

(W) m ° c = exp(-aS), (16) 
where the string tension a is given by the relation 



7T 



<^=^J dz 2 /^). (17) 

The value of the string tension a ~ 0.18 GeV 2 may be determined phenomenologically from the 
slope of the meson Regge trajectories ||. Note that the value T g w 0.13 fm can be extracted 
from the gluelump spectrum, which is calculated in MFC \T2\ using the only nonperturbative 
parameter o. 

From (p~3| ) and (|TJ]) the equation 

* = ^ (18) 



follows, which relates in particular the gluon condensate (^-F^ u (0)F^ u (0)^ = (^F 2 ^ with the 
value of the formfactor D at zero. Although D(0) does not enter formally into the definition 
of cr (|T7D , but assuming the behaviour of -D(z) (0) at small distances (this assumption is 



rigorously speaking false |13| , nevertheless we use it in this paper as an approximation; the 
accurate expressions at small distances can be found in |T^|), one gets a = irD(0)Tg, which 



allows to evaluate the gluon condensate as follows, (^F 2 ) = 18cr/(7r 3 T 2 ) w 0.25 GeV 4 , for 
T g = 0.13 fm. 

The value of the gluon condensate can also be obtained from the ITEP sum rules (see the 



book [T3 1 and references to the original works therein) . The sum rules use the operator expan- 
sion for the currents, and the gluon condensate contribute to it in combination {^G 2 } /Q 4 . 
Herewith the value of the gluon condensate is assumed to be independent on the momen- 
tum. Recently on the basis of the analysis using the sum rules of the experimental data on 
r-lepton decays and charmonium spectrum in the papers |Tj| were obtained the following 



values of gluon condensate: <^G 2 > = (0.006 ± 0.012) GeV 4 and <^G 2 ) = (0.009 ± 0.007) 
GeV 4 correspondingly. 

One can assume that the sum rules operate with the some effective value of the condensate, 
which equals to the correlator \^F® (x e i)F® (0)) at some effective separation x e f ^ 0. Taking 
(^G 2 ) = 0.007 GeV 4 and r/w 0.13 fm, one obtains x ei = T g \n((fF 2 )/(^G 2 )) w 
0.5 fm, i.e. characteristic radius of confinement. Note that when assuming that bilocal 
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correlators do not depend on point separation and -D(x)=const, and integrating in (fP7|) up to 
the characteristic hadronic size = 1 fm, one can estimate the gluon condensate from (|T7|), 
flIB"|) as (^F 2 ^ ps cr/r^ ps 0.007 GeV 4 , which is also consistent with the value obtained in JlfJ. 

The conclusion is that the sum rules results can not be used neither as a confirmation nor 
as a disclaimer of the MFC parameterization fll4"|). 



3 The baryon potential 

According to (Bp, flI2p, ( |T4D the baryon potential reads as 

y B = liml^ / f da^^da^i^V^^x - aO- 

— i 2 3 a a 



a=l 



'\Yil /^(^SU^mi^^-^). (19) 

where dcr^ denotes an integration over the surface 5 a . Note that the nondiagonal part of the 
potential (i.e. second term in ( |TP| ) at a ^ b) differs in the factor —1/2 = —1/{N C — 1) from 
the corresponding quantity in [18], where the error was made by the author (in eqs.(24) 
from and (46) from pfl). 

Since the surfaces of the Wilson loop are oriented along the temporal axis, only correlators 
of the color-electric field contribute to the potential (|19D , and one gets 

V-V nj a) nf / / rfZrfZ' / dtV Uij 4(z ab ), (20) 

where n( a ) is a vector of unit length directed along the line connecting string junction with 
the corresponding quark; z a \, = (ln^ — Z'n^,t), R a is a distance from the string junction to 
the corresponding quark, and 

VuM-S ij D(z) + -^-^^. (21) 

In what follows we neglect the formfactor Di in (0), since its contribution to the potential 
is small. 

In calculation of the baryon potential one should distinguish configurations of quarks, in 
which quarks form 1) triangles with each of the angles less that 27r/3 and 2) triangles having 
one angle greater than 2tt/3. In the former case the minimal surface of the Wilson loop 
consists of three surfaces intersecting with each other at angles 27r/3, as is shown in Fig. 
1. In the latter one the string junction position coincide with the position of quark at large 
angle, and the Wilson loop consists of two surfaces. Then one can set in Eq. fl20D R 3 = 0. 

Consider in advance configuration of the equilateral triangle. In this case Ri — R2 — R3 
I! . and baryon potential reads as 

V B (R)=3V M (R) + V nd (R), (22) 

where Vm is diagonal and V nd nondiagonal terms. From eqs.(0), (pO]) , (21) one gets 



V M {R) = 2D(0) / dz 1 (R-z 1 ) f dtexpf-^ 
Jo Jo \ L g 
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Figure 2: The lattice nonperturbative baryon potential from || (points) for lattice parameter 
(3 = 5.8 and MFC potential V {B) (solid line) with parameters a = 0.22 GeV 2 and T g = 0.12 
fm vs. the minimal length of the string L. The dotted line is a tangent at L = 0.7 fm. 



= ^ jf^ dxxK^x) -T 3 ( 2 -^ 2 g))l (23) 

where a = nD(0)Tg, K\ and are McDonald functions. This potential determines an 
interaction of static quark and antiquark at separation R. Note that it is related with the 
gluon field in meson calculated in MFC [|18j using the connected probe as follows, 



dV M (R) 2a f R/T '- 



dR 



2o- r n/Is 

— / dxxKx(x) = E (R), (24) 

7T Jo 



where E (R) is a value of the confining field acting on the quark. 
The nondiagonal potential, 

VU*, = *T, - ^-f f ^A, f | , (25, 

2-k T g J* cos ip \ 2T g cosip ' 

is positive, rises from zero to the value 2/\/ r 3aT g , and saturates at i? > 0.6 fm. 

The behavior of the potential (^) versus the minimal length of the baryon string L = 
Ri + i?2 + Rz — 3-R is shown in Fig. 2 in comparison with the lattice data from 0. One can 
see that the potential completely describes the lattice results. At large distances L > 1.5 fm 
the potential rises linearly with the slope a. At characteristic hadronic sizes the slope of the 
potential falls, which is shown in Fig. 2 by the tangent at point L = 0.7 fm with the slope 
~ 0.9cr. Note that all lattice points at L < 1.5 fm are cosistent with the tangent. This effect 



is also confirmed by the phenomenology of baryon spectrum |L9 1 . 

To perform the analysis of the lattice results from 0, one should add to the potential ( p2]) 
the perturbative color-coulomb one, 
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Figure 3: The lattice baryon potential in the equilateral triangle with quark separations r 
from § (points) at /3 = 5.8 and the MFC potential + V£*L (solid line) at a s = 0.18, 



a = 0.18 GeV 2 , and T g = 0.12 fm. 



where Cp = 4/3 is the fundamental Casimir operator. The results are shown in Fig. 3. One 
can see the full agreement of our potential with this independent set of lattice data as well. 

In the case of the triangle with different separations R a of quarks from the string junction, 
having angles no larger than 27r/3 



v b {r u r 2 , r 3 ) = v M {R a ) + V ^{Ra, Rb), 



(27) 



a=l 



i<b 



where 



<Pab 



dip 



VSR a 



\- IP n , ^II _ J p-2 



+ 



+4 f' -%-k, ^ 



(28) 



If one of the angles of the triangle is larger than 27r/3, the position of string junction 
coincides with the position of quark at large vertex. Let us denote a an angle supplementary 
to the large one; r 1; r 2 the distances from the string junction to the adjacent quarks. Then 

we get 

V B (r u r 2 , a) = V^n) + V M (r 2 ) + V nd { ri , r 2 , a), (29) 

where 



V nd (ri,r 2 ,a) 



2a ctg a 



71 



+r 2 2 



sin 2a a 
2ixT n 



d ip 



-Ko 



r± sm a 



v sm z (p \T g sm<p 



o sin 2 (a — if) \T g sin(a — tp) 
r 2 sin a x 



+ 



(30) 



Note that Kd^i, r 2 , a) = V nd (r 2 , r x , a); V vd (r 1 ,r 2 , 7r/3) = V nd (ri, r 2 ), V nd (R, R) = 1/3! V nd (R). 
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For the analysis of the dependence of potential on the location of quarks, let us consider 
its behavior in isoceles triangles with the vertex 0<7<<27r/3 and fixed length of the string 
L = Ri + R 2 + -R3 versus the value of 7. At distances R a > 0.5 fm one can use asymptotic 
formulas following from (p3[) , (p8f) and (|30|), 



V M (R)aaR--oT g , 

7T 

V Qd (R a ,R b ) » ^f, V nd (r u r 2 ,a) ^^^aT g . (31) 
3v 3 TT 

One can conclude that at 7 = 0, when positions of two quarks coincide and these quarks 
are in the antitriplet of color 577(3), the baryon potential in isoceles triangle reduces to the 
meson one, 

V L (j = 0) « ah - - aT g . (32) 

7T 

If < 7 < 27r/3, then the string consists of three lengths and 

V L (0 < 7 < 2tt/3) « aL + f-H + aT 5 . (33) 
If the string consists of two lengths, then 

V L ( 7 > 2tt/3) « + f -- + -(tt - 7) ctg (vr - 7) J crT„. (34) 
According to eqs. (|32| ) - (^4|), when 7 icreasing, the potential rapidly falls by the value 

AV 1 = (-- --L ] aT 5 w 150 MeV, (35) 



v/37 



almost does not change in the range < 7 < 27r/3, then rapidly increases at 7 ~ 27r/3 by the 
value 

AV 2 =(i- aT 9 « 55 MeV, (36) 

and slowly increases in the range 27r/3 < 7 < tt by the value 



AVi . K» - V\j) = ^ - ^ j .T 5 « 30 MeV, (37) 

where numerical values are given for a = 0.18 GeV 2 , T g = 0.12 fm. The dependence of 
the baryon potential in isoceles triangle on the angle 7 at L = 1.8 fm is shown in Fig. 4. 
Note that character jumps and drops of the potential are proportional to the quantity cT g , 
i.e. determined by confinement as well as by correlations of the stochastic nonperturbative 
fields, and it would be interesting to obtain them independently in lattice QCD. Baryon 
potential in various configurations of quarks was considered on the lattice in ||. However, 
configurations with the angles in the range 7r/20 < a < n/2 considered in this work did 
not allow to establish this dependence, because in the range of the small angles the size of 
lattice spacing becomes comparable with the quark separation and accurate calculations of 
the potential become difficult. The further lattice investigations of configurations with large 
angle are called for. 
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Figure 4: The baryon potential in isoceles triangle with the length of the string L = 1.8 fm 
versus the vertex 7 for a = 0.18 GeV 2 , T g = 0.12 fm. 



It is also interesting to consider the energy density of the string near the string junction, 
or the quantity 



V L (i) 



l - V M {Ra) -aA+Y, V nd(Ra, Rb). (38) 

\a=l / a<b 

If the vertex in the isoceles triangle with the string length fixed is in the region 2tc/3 < 7 < n, 
the potential v (7) has the same behavior as the V 1 ^), since in this region the dependence 
of the latter on the value of vertex is determined by the only nondiagonal term. This means 
that while the vertex is increasing in this range, the energy density is rising. The drop at 
7 = for V is small, 

AVj. = f ^ - -j=) aT 9 = 0-12aT 9 « 13 MeV. (39) 

At 7 = 27r/3 the potential V L (j) drops but not jumps as V L {j) does, 

AV 2 = - cxT g = -0.13aT g w -15 MeV, (40) 

i.e. the energy density decreases. One can conclude that the energy density near the string 
junction depends on the location of the quarks only slightly. This is in contradiction with the 
field distributions obtained in |17|, |18[ using the connected probe, where author has made an 
error in the basic formulas (11)-(18) from [17] and (27) from [fT8]| . It is not hard to show that 
the cited formulas and corresponding field distributions concern to the two or three static 
meson Wilson loops in the special case when the positions of antiquarks coincide, but not to 
baryons. The problem of the field distributions in baryon with the connected probe has some 
ambiguities and will be considered in subsequent publications. 
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4 Conclusions 



In the present work the nonperturbative static potential in baryon was calculated analytically 
in the framework of MFC for arbitrary locations of quarks. Two important effects related 
with the correlations of the nonperturbative gluon fields are studied: a decrease of the slope of 
the potential at hadronic lengths, and its dependence on the locations of quarks when length 
of the string fixed. The latter effect is found for the first time and shown to be proportional 
to the combinations of parameters crT g . It would be interesting to verify it independently on 
the lattice and extract an accurate value of the correlation length of gluonic fields. 



The errors in preceding works JL7| [[H were corrected for baryon potential and stated for 
the field distributions in baryon with the connected probe. 

The dependence of the potential on the length of the string was shown in the paper to be 
in complete agreement with the precise numerical lattice calculations, which allows to perform 
the accurate analytic studies of the baryon spectra, first of all of the spectrum of nucleons, 
and has important meaning for various applications in nuclear physics. 

The author is grateful to Yu.A.Simonov for numerous useful discussions. I also thank 
T.Takahashi, Ph.de Forcrand and V.I.Shevchenko for the correspondence. Partial support by 
RFBR grants 00-02-17836, 00-15-96786, and INTAS 00-00110, 00-00366 is acknowledged. 
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